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Abstract 

The theoretical status of Higgs boson and supersymmetric particle production at 
hadron colliders is reviewed with particular emphasis on recent results and open 
problems. 


1 Introduction 

Supersymmetry imposes a new symmetry between the fermionic and bosonic degrees of 
freedom |l| . Since supersymmetric theories do not develop quadratic divergences in higher 
orders, they provide a natural solution of the hierarchy problem at the electroweak scale 
p. If supersymmetric grand unihed theories (SUSY-GUT) are considered, the predicted 
value of the Weinberg angle turns out to be in excellent agreement with the present 
measurements at the LEP and SLC experiments 0. Moreover, owing to the large top 
quark mass SUSY-GUTs develop electroweak symmetry breaking at the electroweak scale 
dynamically [Q. Due to these properties the supersymmetric extension of the Standard 
Model exhibits one of the most attractive alternatives beyond the Standard Model. 

The Higgs mechanism is a cornerstone of the Standard Model (SM) and its supersym¬ 
metric extensions |P. Thus, the search for Higgs bosons is one of the most important 
endeavors at future high-energy experiments. The minimal supersymmetric extension of 
the Standard Model (MSSM) requires the introduction of two Higgs doublets in order to 
preserve supersymmetry. There are hve elementary Higgs particles, two GP-even {h,H), 
one GP-odd (A) and two charged ones {H^)- At lowest order all couplings and masses 
of the MSSM Higgs sector are hxed by two independent input parameters, which are 
generally chosen as tan/3 = V 2 /V 1 , the ratio of the two vacuum expectation values 
and the pseudoscalar Higgs-boson mass M 4 . At LO the light scalar Higgs mass Mh has 
to be smaller than the Z-boson mass Mz- Including the one-loop and dominant two- 
loop corrections the upper bound is increased to Mh ^ 135 GeV 0. The negative direct 
searches for the Higgsstrahlung processes > Zh, ZH and the associated production 

e+e“ —> Ah, AH yield lower bounds of rrih^H > 91-0 GeV and rriA > 91.9 GeV. The 
range 0.5 < tan/3 < 2.4 in the MSSM is excluded by the Higgs searches at the LEP2 
experiments [ 0 . 

* Talk presented at SUSY02, DESY Hamburg, 2002. 
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Higgs bosons can be searched for at the upgraded Tevatron, a pp collider with a c.m. 
energy of 2 TeV, and the LHC, a pp collider with a c.m. energy of 14 TeV. At the Tevatron 
the most important processes are Higgs-strahlung qq ^ W + h/H with h/H bb which 
is important for Higgs masses below about 130 GeV, gluon fusion gg ^ h/H ^ W*W 
which is relevant for Higgs masses above 130 GeV, and Higgs radiation off bottom quarks 
qq —> bb(j) which plays an important role for large values of tan f5. With an integrated 
luminosity of 30 fb“^ the Tevatron can probe the entire MSSM parameter space [^. At the 
LHG the most important Higgs production modes are gluon fusion, vector boson fusion 
qq ^ qq + h/H and Higgs radiation off top and bottom quarks. There are several Higgs 
decay modes which enable the discovery of the Higgs bosons 0. 

The novel colored particles, squarks and gluinos, and the weakly interacting gauginos 
can be searched for at the Tevatron and the LHG. Until now the search at the Tevatron 
has set the most stringent bounds on the colored SUSY particle masses. At the 95% 
GL, gluinos have to be heavier than about 180 GeV, while squarks with masses below 
about 180 GeV have been excluded for gluino masses below ~ 300 GeV Stops, the 


scalar superpartners of the top quark, have been excluded in a signihcant part of the 
MSSM parameter space with mass less than about 80-100 GeV by the LEP [|^ and 
Tevatron experiments |OT. Finally charginos with masses below about 100 GeV have 
been excluded by the LEP experiments [^, while the present search at the Tevatron is 


sensitive to chargino masses of about 60-80 GeV with a strong dependence on the specihc 
model . Due to the negative search at LEP2 the lightest neutralino Xi has to be heavier 
than about 45 GeV in the context of SUGRA models 0. In the i?-parity-conserving 
MSSM, supersymmetric particles can only be produced in pairs. All supersymmetric 
particles will decay to the lightest supersymmetric particle (LSP), which is most likely 
to be a neutralino, stable thanks to conserved i?-parity. Thus the hnal signatures for the 
production of supersymmetric particles will mainly be jets, charged leptons and missing 
transverse energy, which is carried away by neutrinos and the invisible neutral LSP. 


2 Higgs boson production 


2.1 Gluon fusion 

The gluon fusion mechanism gg ^ (f) provides the dominant production mechanism of 
Higgs bosons at the LHG in the entire relevant mass range up to about 1 TeV for small 
and moderate values of tan/? in the MSSM [|^. At the Tevatron this process plays a role 


for Higgs masses between about 130 GeV and 190 GeV, if the branching ratio of decays 
into W*W pairs is large enough |^. The gluon fusion process is mediated by heavy quark 
triangle loops and, in the case of supersymmetric theories, by squark loops in addition, if 


the squark masses are smaller than about 400 GeV |T4 


In the past the full two-loop QGD corrections have been determined. They increase 
the production cross sections by 10-90% [|T^, thus leading to a significant change of 
the theoretical predictions. Very recently, the full NNLO calculation has been hnished 

0 


in the heavy top quark limit 


This limit has been demonstrated to approximate 


the full massive K factor at NLO within about 20% for small tan/? in the entire mass 
range up to 1 TeV [0. Thus, a similar situation may be expected at NNLO. The reason 
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for the quality of this approximation is that the QCD corrections to the gluon fusion 
mechanism are dominated by soft and collinear gluon effects, which do not resolve the 
one-loop Higgs coupling to gluons. Fig. ^ shows the resulting iF-factors at the LHC and 
the scale variation of the iF-factor for the SM Higgs boson. The calculation stabilizes 
at NNLO, with remaining scale variations at the 10-15% level. These uncertainties are 
comparable to the experimental errors which can be achieved with 300 fb“^ of data at 
the LHC 1^. The full NNLO results conhrm earlier estimates which were obtained in 
the frame work of soft gluon resummation and soft approximations of the full 


three-loop result within 10-15%. The full soft gluon resummation has been performed in 
Ref. 0. The resummation effects enhance the NNLO result further by about 10% thus 


signaling a perturbative stabilization of the theoretical prediction for the gluon-fusion 
cross section. 



Figure 1: Scale dependence of the K-factor at the LHC. Lower curves for each pair are 
for = 2 Mh, /iF = Mh/ 2, upper curves are for = Mh/ 2, p,p = 2 Mh. The K-factor 
is computed with respect to the LO cross section at = jap = Mp. From Ref. 0/. 


In supersymmetric theories the gluon fusion cross sections for the heavy Higgs, H, and, 
for small Ma, also for the light Higgs, h, are signihcantly affected by bottom quark loops 
for tan /? > 3 so that the heavy top quark limit is not applicable in general. This can be 
clearly seen in the NLO results, which show a decrease of the K factor down to about 1.1 
for large tan (3 H- This decrease originates from an interplay between the large positive 
soft/collinear gluon effects and large negative double logarithms of the ratio between the 
Higgs and bottom masses. In addition, the shape of the pp distribution of the Higgs 
boson may be altered; if the bottom loop is dominant, the pp spectrum becomes softer 
than in the case of top-loop dominance. These effects lead to some model dependence of 
predicted cross sections. 


2.2 tt(f) production 

SM Higgs boson production in association with tf pairs plays a signihcant role at the 
LHC for Higgs masses below about 130 GeV, since this production mechanism makes 
the observation of if 66 possible |^. The decay iL ^ 77 is potentially visible in 
this channel at high integrated luminosity. For Higgs masses above about 130 GeV, the 
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decay H W*W can be observed [^. ttH production could conceivably be used to 


determine the top Yukawa coupling directly from the cross section. NLO QCD corrections 
have become available. They decrease the cross section at the Tevatron by about 20% 
1^, while they increase the signal rate at the LHC by about 20-40% H, see Fig. |. 


The scale dependence of the production cross section is signihcantly reduced, to a level 
of about 10-15%, which can be considered as an estimate of the theoretical uncertainty. 
The transverse momentum and rapidity distributions at NLO can be approximated by a 
rescaling of the LO distributions with a constant K factor within 10-15% [^. Thus, the 
signal rate is under proper theoretical control now. In the MSSM, tth production with 
h —> 77 , bb is important at the LHC in the decoupling regime, where the light scalar h 
behaves as the SM Higgs boson ill- Thus, the SM results can also be used for tth 
production in this regime. 



M„ [GeV] 


Mh [GeV] 


Figure 2: The cross section for pp/pp —> ttH + X at the Tevatron and LHC in LO 
and NLO approximation, with the renormalization and factorization scales set to p = 
mt + Mh/2. 


2.3 bb(j) production 


In supersymmetric theories bb(j) production becomes the dominant Higgs boson production 
mechanism for large values of tan (3 [|T^ , where the bottom Yukawa coupling is strongly 
enhanced. In contrast to ticf production, however, this process develops potentially large 
logarithms, logm^/m^, in the high-energy limit due to the smallness of the bottom quark 
mass, which are related to the development of b densities in the initial state. They 
can be resummed by evolving the b densities according to the DGLAP-equations and 

B 


introducing them in the production process 
5-initiated processes bb ^ H and b g 


^b^H 


The NLO QCD corrections to the 
] are known to be moderate. 
The resummation increases the cross section by a factor of about 5 at the Tevatron 
and about 2-3 at the LHC and thus plays a significant phenomenological role. The 
introduction of conventional b densities, however, requires an approximation of the hard 
process kinematics, i.e. the initial and hnal b quarks are assumed to be massless and travel 
predominantly in forward and backward direction. These approximations can be tested 
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in the full gg —> bbcp process. 

We have to investigate if the energy of the Tevatron and LHC is sufficiently large to 
develop the factorization of bottom densities, i.e. that the transverse mass distribution of 
the (anti)bottom quark can be factorized as a convolution 

dcr I (a s T-, - 1 -1 \ 

—— = —Pqg ^g^g^abg} + non-smgular terms (I) 

alVlTb Jypb I J MTb=mb^0 

where Mxb = \J + Prb denotes the transverse mass of the (anti)bottom quark, Pgg the 
corresponding DGLAP splitting kernel, g the gluon density of the (anti)proton and abg 
the partonic cross section for ^b^g H. This factorization requires that the transverse 
mass distribution is dominated by the hrst term, i.e. da/dMrpb oc l/Mrpb-, for transverse 
masses up to the factorization scale of the (anti)bottom density. The transverse mass 

-3 




[GeV] Mt, [GeV] 

Figure 3: Transverse mass distributions of the bottom quark in bbH production at the 
Tevatron and the LHC. We have adopted CTEQ5M1 parton densities and a bottom mass 
of mb = 4.62 GeV. The solid lines show the full LO result from qq^gg bbH and the 
dashed lines the factorized collinear part of Eq. W which is absorbed in the bottom parton 
density. The upper curves are multiplied with the factor MTb of the asymptotic behavior, 
which is required by factorizing bottom densities. 

distributions at the Tevatron and LHC are shown in Fig. The solid curves show 
the full distributions of the qq, gg bbcf processes, while the dashed lines exhibit the 
factorized collinear part of Eq. o which is absorbed in the bottom density. For a proper 
factorization, these pairs of curves have to coincide approximately up to transverse masses 
of the order of the factorization scale, which is usually chosen to be pip = 0{mH)- It is 
clearly visible that there are sizeable differences between the full result and the factorized 
part, which originate from sizeable bottom mass and phase space effects, that are not 
accounted for by an active bottom parton density. Moreover, the full result falls quickly 
below the approximate factorized part for transverse masses of the order of which 
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is much smaller than the usual factorization scale used for the bottom densities. We 
conclude from these plots that bbcj) production at the Tevatron and LHC develops sizeable 
bottom mass and kinematical phase space effects, so that the use of bottom densities in 
the process 66 —> 0 may lead to an overestimate of the correct theoretical result due to too 
crude approximations in the kinematics of the hard process |^. The full NLO calculation 
of the gg —>■ bh(f> will yield much more insight into this problem, since the large logarithms 
related to the evolution of bottom densities have to appear in the NLO corrections, if the 
picture of active bottom quarks in the proton is correct. 


3 SUSY particle production 


3.1 Production of squarks and gluinos 

Squarks and gluinos can be produced via pp,pp —> qq, qq, qg, gg at hadron colliders. The 
determination of the full SUSY-QCD corrections has been performed for the upgraded 
Tevatron and the LHC. For the natural renormalization/factorization scale choice Q = m, 
where m denotes the average mass of the hnal-state SUSY particles, the SUSY QCD 
corrections are large and positive, increasing the total cross sections by 10-90% a. 
The inclusion of the NLO corrections reduces the LO scale dependence by a factor 3-4 
and reaches a typical level of ~ 10-15% which serves as an estimate of the remaining 
theoretical uncertainty. Moreover, the dependence on different sets of parton densities is 
rather weak and leads to an additional uncertainty of ~ 10-15%. In order to quantify the 
effect of the NLO corrections on the search for squarks and gluinos at hadron colliders, 
the SUSY particle masses corresponding to several hxed values of the production cross 
sections have been extracted. These masses are increased by 10-30 GeV at the Tevatron 
and 10-50 GeV at the LHG, thus enhancing the present and future bounds on the squark 
and gluino masses signihcantly. Finally, the QGD-corrected transverse-momentum and 
rapidity distributions for all different processes have been evaluated. The modihcation 
of the normalized distributions in NLO compared to LO is less than about 15% for the 
transverse-momentum distributions and much less for the rapidity distributions. Thus it 
is a sufficient approximation to rescale the LO distributions uniformly by the K factors 
of the total cross sections . 


3.2 Stop pair production 

At LO only pairs of U or pairs of 62 can be produced at hadron colliders. QGD-initiated 
mixed tiO pair production is only possible at NLO and beyond. However, mixed stop pair 
production is completely suppressed by several orders of magnitude and can thus safely 
be neglected |^. The evaluation of the SUSY-QGD corrections proceeds along the same 
lines as in the case of squarks and gluinos. They increase the total cross sections by up 
to about 40% [0. As in the squark/gluino case the scale dependence is strongly reduced 


and yields an estimate of about 10-15% of the remaining theoretical uncertainty at NLO. 
At NLO the virtual corrections depend on the stop mixing angle, the squark, gluino and 
stop masses of the other type. However, it turns out that these dependences are very 
weak and can safely be neglected. 
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3.3 Chargino and neutralino production 


The production cross sections of charginos and neutralinos depend on several MSSM 
parameters, i.e. Mi,M 2 ,/i and tan/? at LO [^]. The cross sections are sizeable for 
chargino/neutralino masses below about 100 GeV at the upgraded Tevatron and less 
than about 200 GeV at the LHG. Due to the strong dependence on the MSSM parame¬ 
ters the extracted bounds on the chargino and neutralino masses depend on the specific 
region in the MSSM parameter space ig. The outline of the determination of the SUSY- 
QGD corrections is analogous to the previous cases of squarks, gluinos and stops. The 
corrections enhance the production cross sections of charginos and neutralinos by about 
10-40% 1^. The scale dependence is reduced to about 10% at NLO, which signalizes 
a signihcant stabilization of the theoretical prediction for the production cross sections. 
The dependence of the chargino/neutralino production cross sections on the specific set 

m- 


of parton densities ranges at about 10-15% 


3.4 Associated production of gluinos and gauginos 


The cross sections of the associated gluino-gaugino production are sizeable for the lightest 
chargino/neutralino states at the upgraded Tevatron and the LHG if the gluino mass is 
less than about 400-500 GeV [^. The determination of the SUSY-QGD corrections is 
analogous to the previous cases of squarks, gluinos, stops and gauginos. The production 
cross sections are decreased by up to 3% at the Tevatron and increased by about 10-20% 
at the LHG due to these corrections as can be inferred from Fig. B The scale de- 




m~ [GeV] m| [GeV] 

Figure 4: K factor of the cross sections for gluinos produced in association with the gaug¬ 
inos X 2 xt O'i the upgraded Tevatron (left) and the LHC (right). Parton densi¬ 
ties: CTEQfL (LO) and CTEQfM (NLO) with the renormalization/factorization scale 
Q = {mg + m^)/2. 


pendence is shown in Fig. ^ at LO and NLO. It is clearly visible that it is reduced to about 
10-15% at NLO which signalizes a significant stabilization of the theoretical prediction 

^ These results disagree with the corresponding figures in Ref. Q. The discrepancies have been 
resolved meanwhile and the corrected calculations in the Erratum to Ref. |Q are in agreement with the 
results that we have presented to this conference. We are grateful to the authors of Ref. Q for their 
cooperation. 
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Q/m Q/m 

Figure 5: Scale dependence of the cross sections of gluinos produced in association with 
X 2 O't the upgraded Tevatron (left) and the LHC (right). Parton densities: CTEQfL (LO) 
and CTEQ4M (NLO) with the renormalization/factorization scale Q varied in units of 
the average mass m = {mg + m/)/2. 


for the production cross sections. The dependence of the gluino-gaugino production cross 
sections on the specihc set of parton densities ranges at about 10-15% [^. 


4 Conclusions 


Considerable progress has been made recently in improving QCD calculations for Higgs 
and supersymmetric signal cross sections at hadron colliders. Most (N)NLO [SUSY-]QCD 
corrections to all relevant production processes at hadron colliders are known, i.e. the 
theoretical status of novel particle production at the Tevatron and LHC is nearly complete. 
Large corrections to many processes have been found which underlines the importance of 
including them in realistic experimental analyses. After inclusion of the NLO corrections 
the residual theoretical uncertainties are reduced to a level of 10-15%. There are several 
Fortran programs which include most of the higher order corrections [^ . 
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